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ABSTRACT 

In this paper we study the "standardized candle method" using a sample of 
37 nearby (redshift z < 0.06) Type II plateau supernovae having BVRI photom- 
etry and optical spectroscopy. An analytic procedure is implemented to fit light 
curves, color curves, and velocity curves. We find that the V — I color toward 
the end of the plateau can be used to estimate the host-galaxy reddening with a 
precision of cr(Ay) = 0.2 mag. The correlation between plateau luminosity and 
expansion velocity previously reported in the literature is recovered. Using this 
relation and assuming a standard reddening law (Ry = 3.1), we obtain Hubble 
diagrams in the BV I bands with dispersions of ~ 0.4 mag. Allowing Ry to vary 
and minimizing the spread in the Hubble diagrams, we obtain a dispersion range 
of 0.25-0.30 mag, which implies that these objects can deliver relative distances 
with precisions of 12-14%. The resulting best-fit value of Ry is 1.4 ± 0.1. 

Subject headings: stars: supernovae: general — cosmology: distance scale - 
galaxies: distances and redshifts 
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INTRODUCTION 



Core-collapse SNe (C CSNe) are in general clos ely associated with star-forming regions 



in late-type galaxies (e.g., lAnderson fc Jamesl 120081 ) . Therefore, they have been attributed 



to stars born with initial mass > 8 M that undergo the collapse of their iron cor es after a 



Burrows 



few m illion years of evolution, and the subsequent ejection of their envelopes (e.g., 
20001 ) . According to theoretical stu dies, these SNe leave a compact object as a remnant , 
either a neutron star or a black hole (IBaade fc Zwickylll934l ; lArnettlll996l ; iJanka et al.ll2007l ). 



Among CCSNe, we can observ ationally distingu ish those with prominent hydrogen lines 
in their spectra (dubbed Type II; iMinkowskil Il94ll ). those with no H but strong He lines 
(Type lb), and those lacking H or He lines (Type Ic); see (jFilippenkol 119971 ) and references 
therein. Although all of these objects are thought to share largely the same explosion 
mechanism, their different observational properties are explained in terms of how much of 
their H-rich and He-rich envelopes were retained prior to explosion. When the star explodes 
with a significant fraction of its initial H-rich envelope intact, it should display a H-rich 
spectrum and a light curve characterized by an optically thick phase of ~ 100 days of nearly 
constant luminosity called a "plateau," which comes to an end at the point when all of the 
hydrogen in the envelope has recombined. This phase is followed by a sudden drop of 2-3 
mag and thereafter by an exponential decay in lumin osity caused by the radioactive decay 
of 56 Co into 56 Fe (e.g., lNadvozhinll20~0~3l lutrc^EA Nearly 50% of all CCSNe belong t o 
this class of "Type II plateau" SNe (SNe IIP; ICappellaro et alill999l : iBotticella et~aDl2008h . 



A decade ago, the standardized lumino sities of Type la SNe (e.g.. |Phillipall993uHamuy et al. 



19961 ; iRiess et al.lll996t iPhillips et al.lll999l ) led to the construction of Hubble diagrams (HDs) 
at redshifts z = 0-0.1 with dispersions of only ~ 0.15 mag, which opened up the opportu- 
nity to measure t he ex pected deceleration of t he Universe. At the end of the millennium, 
Riess et al.l (119981 ) and iPerlmutter et al.l (Il999l ) reported independent observations of high- 
redshift SNe (z = 0.1-0.5) which permitted the measurement of the expansion history of 
the Universe over ~ 5 Gyr of lookback time. Contrary to expectations, these observations 
revealed that the Universe is currently described by an accelerated expansion. The discovery 
of co smic acceleration, which was subsequently confirmed by additional studies of SNe l a 
(e.g., lAstier et alJbood : IWood-Vasev et al.ll2007t iHicken et allbood : iFreedman et~aD 120091 ) . 



implies the possible existence of a cosmological constant, a concept initially introduced by Al- 
bert Einstein at the beginning of the 20th century. Its origin is still a mystery, and there are 
man y other possible can didates for the repulsive "dark energy" that permeates the Universe 



sec 



Frieman et al. 2008 for a review 



The acceleration of the Universe has now been indirectly confirmed by o ther indepen- 
dent experiments such as the Wilkinson Microwave Anisotropy Probe (WMAP; ISpergel et al. 
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Boughn fc Crittenden 



200 71; iBennett et al.ll2003f) and baryon acoustic oscillations (BAO; lBlake Glazebrookll2003 
Seo fe Eisensteinll2003), as well as more directly with the integrated Sachs- Wolfe effect (e.g 



2004 iFosalba et al.ll2003l ) and X-ray observations of clusters of galax- 



ies ( lAUen et al.l 120081 ). but it is nevertheless important to obtain additional independent 
confirm ation of the SN l a resul ts. Although SNe IIP are not as luminous and uniform as 
SNe la, lHamuy fc Pintol (120021 ) showed that the luminosities of the former can be stan- 
dardized to levels of 0.4 and 0.3 mag in the V and / bands, respectively. This finding 
converted these objects into precise distance indicators and potentially into useful tools to 
measure cosmological parameters. Moreover, SNe IIP would be more abundant at high red- 
shift than SNe la, since the star-formation rate increases with redshift until z ~ 2.5 (see 
Hopkins &: Beacomll2006l . and references therein). 



Type IIP supernovae provide tw o routes to distance det erminations. First, the "ex- 
panding photosphere method" (EPM; iKirshner fc Kwanlll974l ) is a tec hnique based on the 



oretical models that is independent of the extragalactic distance scale (jEastman et al.lll996 



Dessart fc Hillierl 120051 ) . In EPM the spectral energy distribution of the SN is approximated 
as a black body, and atmosphere models permit one to correct the dilution of the photospheric 
flux relative to a Planck function d ue to electron scatterin g. From a detailed analysis of the 
properties of the SN atmospheres, Eastman et al.l (119961 ) showed that the most important 
variable in determining the dilution factors is the effective temperature. This result permits 
one to determine the photospheric angular radius of the SN from a measurement of the color 
temperature, and infer the distance from the expansion velocity measured from spectral lines, 
without the need to craft specific models for each SN. This method can achiev e dispersions of 



0.3 mag in the HP, which translates in to a 14% precision in relative distances (jSchmidt et al. 



1994; 


Hamuv 


2001; 


Jones et al. 


2009) 



ing corrections and color temperatures from tailored models for eac h observation, reyealing 
an accuracy of ~ 1 0% in distance determination of individual SNe (IDessart fc Hillierl 12006 



Dessart et al J 120081 ) . 



Seco nd, the "standardized c andle method" (SCM) is an empirical technique initially pro- 
posed by lHamuy &; Pintol (120021 ) based on an observational correlation between the absolute 
magnitude of the SN and the expansion velocity of the photosphere, the luminosity/expansion- 
velocity (LEV) relation. This correlation shows that SNe IIP with greater luminosities have 
higher expansion velocities, which permits one to correct for the large (~ 4 mag) luminosity 
differences displayed by these objects to levels of only ~ 0.3 mag. So far, the SCM has 
been applied to 24 low - redshift (z < 0.05) SNe (lHamuyll2003bl 120051 ) and more recently by 
Poznanski etall f|2009h to 17 low-redshift (z < 0.03) SNe, as well as bv lNugent etall (bood ) 
to 5 high-redshift (z < 0.29) SNe. The latter work was the first attempt to derive cosmo- 
logical parameters from SNe IIP and demonstrated the enormous potential of this class of 
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objects as cosmological probes. 



Several major sky-survey programs now in progress or to be deployed in the upcom - 



ing years, such as the Palomar Transient Factory (PTF; Rau et al. 2009 ; Law et al. 2009h 



the P anoramic Survey Telescope and Rapid Respons e System (Pan-STA RRS; iHodapp et al. 



20041 ). the Large Synoptic Survey Telescope (LSST; iTyson et al.l 120031 ). the Visible and In- 
frared Survey Telesco pe (VISTA; Emerson et al.l 120041 ) , the VLT Surve y Telescope (VST; 
Capaccioli et al.l 120031). th e Dark Energy Survey (DES; ICastanderl 120071 ) . and the SkyMap- 



per (IGranlund et al. 



20061 ). offer the promise of discovering SNe by the thousands. Thus, we 



will be inevitably confronted by enormous amounts of data on SNe IIP which will contain 
potentially valuable cosmological information. The challenge is to find an implementation of 
SCM that does not require spectroscopic data since these surveys will be purely photometric. 

In spite of the great potential shown by the SCM, it still suffers from a variety of problems 
which need to be addressed: (1) the lack of a well-defined maximum in the light curves has 
prevented us from defining the phase of each event; (2) each SN shows a different color 
evolution, which has compromised the use of the photometric data for t he determination 
of hos t-galaxy extinction; (3) the assumption for dereddening employed by lHamuy fc Pinto 
( 120021 ). namely that all SNe reach the same asymptotic color toward the end of the plateau 
phase, has often led to negative extinction corrections, so it still needs to be further tested; 
and (4) the small sample size used so far, especially the scarcity of SNe in the Hubble flow, 
has prevented a proper determination of the intrinsic precision of the method. 

With these obstacles in mind, we have developed a robust mathematical procedure to 
fit the light curves, color curves, and velocity curves, in order to obtain a more accurate 
determination of the relevant parameters required by the SCM (magnitu des, colors, and 
ejecta velocities), as described in detail in the thesis work of lOlivaresI (120081 ). Here we apply 
this technique to 37 SNe IIP in order to construct a HD, evaluate the accuracy of the SCM, 
and obtain an independent determination of the Hubble c onstant. Since ou r sample has 
several objects in common with the recent EPM analysis of I Jones et al.l (120091 ). we perform 
a comparison between SCM and EPM. 

We organize this paper as follows. In § [2] we describe all of our observational material 
including telescopes, instruments, and surveys involved. The analysis, methodology, and 
procedures, such as the Aq, K, and Ah os t corrections, are explained in detail in § [3j The 
dereddening analysis, the HD, the value of the Hubble constant, and the comparison between 
SCM and EPM distances are addressed in §|H The final remarks in §0 provide a discussion 
about variations of the reddening law for SNe IIP and the possibility of a nonstandard 
extinction law in the SN host galaxies. We present our conclusions in § [61 
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2. OBSERVATIONAL MATERIAL 

This work makes use of data obtained in the course of four systematic SN follow-up 
programs carried out during the years 1986-2003: (1) the Cerro Tololo SN program (1986- 
1996); (2) the Calan/Tololo SN program (CT, 1990-1993); (3) the Optical and Infrared 
Supernova Survey (SOIRS, 1999-2000); and (4) the Carnegie Type II Supernova Program 
(CATS, 2002-2003). As a result of these efforts, photometry and spectroscopy (some infrared 
but mostly optical) was obtained for nearly 100 SNe of all types, 51 of which belong to the 
Type II class. Of these 51 SNe II, a subset of 33 objects meets the requirements for this 
study. Next we describe in general terms the data acquisition and reduction procedures. For 
more details the reader can refer to M. Hamuy et al. (2010, in preparation). 



2.1. Photometric Data 



The photometry was acquired with telescopes from Cerro Tololo Inter-American Ob- 
servatory (CTIO), Las Campanas Observatory (LCO), the European Southern Observatory 
(ESO) in La Silla, and Steward Observatory (SO). Many different telescopes and intruments 
were used to generate this dataset as shown in Table [TJ In all cases w e employed CCD de 



Cousins 


1971; 


Hamuv et al. 


2001) 



The images were processed with IRAP^I through bias subtraction and flat-fielding. All of 
them were further processed through the step of galaxy subtraction using template images 
of the host galaxies. Photometric sequences w ere established around each SN based o n 
observations of Landolt and Hamuy standards (ILandoltl Il992l ; lHamuy et al.l Il992l . Il994l ). 
The photometry of all SNe was performed differentially with respect to the local sequence 
on the galaxy-subtracted images. The transformation of instrumental magnitudes to the 
standard system was done by taking into account a linear color term and a zero point. 
Although this procedure partially removes the instrument-to-instrument differences in the 
SN magnitudes, it should be kept in mind that significant syste matic discrepancies can still 



remain owing to the nonstellar nature of the SN spectrum (e.g., lHamuy et al.lll990l ). 



Complement ary photometry for SN 2003gd obtained by IVan Dyk et al.l ( 120031 ) and 
Hendry et al.l ( 120051 ) was also incorporated in our analysis. 



9 IRAF, the Image Reduction and Analysis Facility, is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., 
under cooperative agreement with the National Science Foundation (NSF); see http : //iraf . noao . edu] 
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2.2. Spectroscopic Data 



The spectroscopic data were also obtained with a great variety of instruments and tele- 
scopes, as shown in Table [TJ The observations consisted of the SN observation immediately 
followed by an emission-line comparison l amp taken at t he same pos ition in the sky, and 2-3 
flux standards per night from the list of lHamuy et al.l (119921 Il994l ). We always used CCD 
detectors, in combination with different gratings/grisms and blocking filters. The reductions 
were performed with IRAF and consisted of bias subtraction, flat-fielding, one-dimensional 
(ID) spectrum extraction and sky subtraction, wavelength calibration, and flux calibration. 
No attempts were made to remove the telluric lines. 

Our spectroscopic database was compl emented with spectra obtained by the CfA SN 
program with the FAST CCD spectrograph (IFabricant et al.lll998l ) at the 1.5-m telescope of 
the Whipple Observatory (see Matheson et all 20081 for details on the data reduction), and by 
the University of Cali fornia, Berkeley SN program with the Kast double CCD spectrograph 
( IMiller &: Stondll993f ) on the Shane 3-m telescope at Lick Observatory. The above standard 
procedures, i ncluding removal of telluric absorption lines, were employed to reduce the Lick 



spectra (e.g.. lFoley et al.ll2003j ) . 



2.3. Subsample Used for this Work 

Of the 51 SNe II observed in the course of these four surveys, a subset of 33 objects 
comply with the requirements of (1) being a member of the plateau subclass; (2) having 
sufficient spectroscopic temporal coverage; and (3) having light curves with good temporal 
coverage. To this sample we added four SNe from the literature: SN 1999gi, SN 2004dj, 
SN 2004et, and SN 2005cs. Table M lists our final sample of 37 SNe IIP. For each SN 
this table includes the name of the host galaxy, the equatorial coordinates of the SN, the 



helio centric redshift and its source, the reddening due to our own Galaxy ( jSchlegel et al. 



19981 ). and the survey during which the SN was observed. 



3. METHODOLOGY AND PROCEDURES 



In the thesis work of lOlivared (120081 ). we describe in detail the techniques for correcting 
the photometry for Galactic absorption (Aq), /^-correction terms, and host-galaxy dust 
absorption (Ah os t) — the so-called U AKA corrections". In brief, we employ a library of 196 
spectra of SNe II (including other SN II types than just SNe IIP) to synthesize iT-terms and 
absorption coefficients given the specific redshift and Ay value of each SN. Synthetic B — V 
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and V — I colors are used as proxies for the spectral energy distribution in order to account 
for the evolution of each SN. 



3.1. Fits to Light, Color, and Velocity Curves 



In the first incarnation of the SCM, lHamuy fc Pinto! ( 120021 ) measured all of the relevant 
quantities (magnitudes, colors, and velocities) at fixed epochs with respect to the time of 
explosion. In most cases, however, it proves hard to constrain this time, thus hampering the 
task to compare data obtained for different SNe. It would be ideal to have a conspicuous 
feature in the light curves; unfortunately, unlike other SN types, SNe IIP generally do not 
show a well-defined maximum during the plateau phase. 

One way around this is to use the end of the plateau as an estimate of the time origin 
for each event. In practice it is not easy to measure this time owing to the generally coarser 
sampling of the light curves at this phase. Thus, our first aim is to implement a robust 
light-curve fitting procedure in order to obtain a reliable time origin to be used as a uniform 
reference epoch to measure magnitudes, colors, and expansion velocities. In the remainder 
of this section, we proceed to implement fitting methods to measure reliable colors and 
expansion velocities. 



3.1.1. BVRI Light- Curve Fits 

Figure [1] shows the BVRI light curves of SN 1999em. As can be seen, there are three 
distinguishable phases in the light curves. 



The plateau phase, in which the SN shows an almost constant luminosity during the first 
~ 100 days of its evolution. This phase corresponds to the optically thick period in 
which a hydrogen recombination w ave recedes in mass coordinate, gradually releasing 
the internal energy of the star (see iKirshner fc Kwanl Il974j ; iNadyozhinl l2003t lUtrobin 



20071 ; iBersten fc Hamuyl 120091 and references therein) . 



The linear or radioactive tail, a linear decay in magnitude (exponential in flux) starting 
about 130 days after the explosion. This phase corres ponds to the optically thi n period 
powered by the 56 Co —> 56 Fe radioactive decay (e.g., IWeaver fc Woosleylll980l ). 



— A transition phase of ~ 30 days between the plateau and linear phases. 



-8 - 



Both the plateau and linear phases are trivial to model if taken separately, but the 
abrupt transition makes the fitting task much more challenging, especially with coarsely 
sampled light curves. After considering several options we ended up using the arithmetic 
sum of the three functions shown in Figure [T] as follows. 



(1) A Fermi-Dirac function (dotted line in Figured]), which provides a very good description 

of the transition between the plateau and radioactive phases: 

/fdW = i + e (*-°p T )M>' (1) 

where a represents the height of the step in units of magnitude, t PT corresponds to 
the middle of the transition phase and is a natural candidate to define the origin of 
the time axis, and wo quantifies the width of the transition phase. At t — tpx — 3u>o 
the height of the step has been reduced by 4.7%, and it decreases down to 95.3% at 
t = tpx + 3wo- 

(2) A straight line (long-dashed line in Figure [I]) which accounts for the slope due to the 

radiactive decay: 

l{t) =po{t- t PT ) + m , (2) 

where po corresponds to the slope of radioactive tail and an approximate slope for the 
plateau in units of mag per day, and m corresponds to the zero point in magnitude 
at t = tpx- 

(3) A Gaussian function (short-dashed line in Figured]), which serves mainly for fitting the 

B light-curve bump and the I light-curve curvature during the plateau phase: 

g(t) = -Pe~^)\ (3) 

where P is the height of the Gaussian peak in units of magnitudes, Q is the center of the 
Gaussian function in days, and R is the width of the Gaussian function. The Gaussian 
function is also useful for reproducing the small-scale features that can appear in the 
K-band plateau. 



The resulting analytic function we use to model the light curves is the sum of the three 
functions detailed above, 

Fit) = f FB (t) + l(t)+g(t) 

-a Q 



X -)_ e (t-tp T )/wo 



+ Po(t- £ PT ) + m - P e~y — ) , (4) 
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which has eight free parameters. This function is fitted to the individual light curves in 
unison (not by parts) using a \ 2 mi nimizing procedure. In order to find the minimum % 2 , we 



use the Downhill Simplex Method ( INelder fc Meadl 119651 ). Despite being not very efficient 



in terms of the number of iterations required, it provides robust solutions. Examples of the 
analytic fits are shown as solid lines in Figure [TJ Although we cannot model most of the 
small-scale features in the plateau, the fitting does a very good job modelling the transition. 
Furthermore, the analytic function gives us important parameters that characterize the light- 
curve shape, particularly tpx which provides a time origin. In the plots in Figure [TJ the time 
axis is chosen to coincide with the value of this parameter obtained from the V light curve. 
In Table [3] we gather the eight parameters and the corresponding reduced \ 2 (xl) °f the 
V light curve fitting for the whole SN sample. 

A critical quantity in the analysis that follows is tpp and its uncertainty, both of which 
determine the uncertainties in all of the relevant SCM quantities. For most cases, when 
the tail phase has been observed, the \ 2 minimizing routine has no difficulties finding tpx 
and delivers a credible error. On the other hand, when the light curve does not have any 
late-time data points, the routine underestimates cr(tpx), in which case we need to provide 
a more realistic estimate of this uncertainty. Only for a few cases did we have to arbitrarily 
set tpx and its error based on the light-curve sampling (see Column 3 in Table [3]). 

Regardless of the sampling of the light curves, we assign a minimum of cr(tpT) = 2 days. 
Although these uncertainties in ipx are very conservative and arbitrary, they do not translate 
into high or unrealistic errors in quantities such as magnitudes, colors, or expansion velocities, 
as we will see later. 



3.1.2. Color-Curve Fits 



Figure |2] shows the B — V, V — R, and V — I colors of three prototypical SNe IIP 
corrected for Aq and i^-terms as described by lOlivaresI ( 120081 . Chapter 3.1). The time origin 
in the abscissa corresponds to tpp, the middle of the transition phase. In each case we 
employ all the data points between days —100 and —10 to fit a Legendre polynomial shown 
with solid lines in Figure |2j The degree of the polynomial was chosen on a case-by-case 
basis and varied between third and sixth order. It is evident that, during the plateau phase, 
the photosphere gets redder with time owing to the decrease of the surface temperature as 
the SN expands and to a redistribution of flux from blue to red wavelengths due to line 
blanketing. According to the theory, the photospheric temperature should approach and 
never get below th e temperature of hydrog en recombination around 5000 K. Based on this 
physical argument, lHamuy fc Pinto! (120021 ) argued that all SNe IIP should reach the same 
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intrinsic colors toward the end of the plateau phase; hence, they proposed that the color 
excesses measured at this phase coul d be attributed to dust r eddening in the SN host galaxy 
and be exploited to measure A^ ost . iHamuy fe Pinto! (120021 ) performed their analysis with 
a simple naked-eye estimate of the asymptotic colors. Here we improve significantly this 
situation through the formal color-curve fitting procedure described above. 

Armed with the polynomial fits, we proceeded to interpolate colors on a continuous grid 
with one-day spacing between days —80 and —10 for analyzing colors at multiple epochs 
(see § 13. 2p . In a handful of cases, the data did not encompass the whole grid and we had to 
extrapolate colors, but never by more than 3 days from the nearest data point. 



3.1.3. Fe II Based Expansion- Velocity Curves 



The third ingredient for SCM is the velocity of the SN ejecta. Different spectral lines 
yield different expansion velocities depending on the region from which the photon is emit- 
ted. The Fe II A5169 li ne is thought to be a good proxy for SN IIP photospheric velocity 
( IDessart fc Hillierl 120051 ) and has been previously em ployed for SCM. Here we use the Fe II 
expansion velocities measured by I Jones et al.l (120091 ) from the minimum of the P-Cygni line 
profile after correcting the spectra for the heliocentric redshifts of the host galaxies. The un- 
certainties of these measurements were estimated to be ~ 85 km s _1 from the quality of the 
wavelength calibration and the signal-to- noise ratio (S/N) of the absorption feature of the 
P-Cygni line profile. Although these measurements do not consider relativistic corrections, 
the latter prove negligible (< 40 km s _1 ) in the range considered here (2000-5000 km s _1 ). 

Figure E] shows Fe II velocities as a function of SN phase, f or four different SNe sel ected 
for their wide range of sampling characteristics. As shown by IHamuy fe Pintol (120021 ). the 
Fe II A5169 expansion velocity curve during the plateau phase can be properly modeled with 
a power law, which we choose to be of the form 



v FeIl (t) = A(t- t ) a , (5) 

where A, to, and a are free fitting parameters without obvious physical meaning. 

In general, we fit for all three parameters, but when only two velocity measurements are 
available (e.g., SN 1992af in the top-left panel of Figure |3]) we fix the a exponent to —0.5, 
which corresponds to a typical value for our sample. As shown with solid lines in Figure |3j the 
fits are quite satisfactory. Here we choose to restrict the power-law fits to the plateau phase, 
since the power-law behavior is not observed for expansion velocities beyond the transition 



- 11 - 



phase. We use the same one-day continuous grid as for the color curves in order to interpolate 
velocities at different epochs between —80 and —10 days. Given the good quality of the fits 
and the shallow slope at late epochs, we allow extrapolations of up to 15 days past the 
nearest data point (see SN 1999cr in Figure [3]). At the early-time boundary, we reduce the 
extrapolations to 10 days prior to the first point, because the power law becomes steeper 
at early times (e.g., see SN 2003gd in Figure |3]). When only two velocity measurements are 
available we reduce the extrapolations by 5 days. Since at the end of the plateau phase the 
SN atmosphere starts to become nebular, we discarded data points beyond the transition in 
the velocity fits. We do not impose any constraints on the first data point. 

3.2. Host-Galaxy Extinction Determination 

To determine the host-galaxy extinction, we assume that all SNe IIP should evolve 
from a hot initial stage to one of uniform photospheric temperature, owing to the hydrogen- 
recombination nature of their photospheres. 

Although simple in theory, there are a few practical difficulties when comparing the color 
evolution. First, it is not always possible to constrain the time of explosion and accurately 
compare color curves of different objects. One way around this is to use the transition time, 
tpT, defined in § 13.1.11 The second problem is illustrated in Figure |5J the color-curve shapes 
can vary significantly from one SN to another, preventing one from measuring a single color 
offset between two SNe. To circumvent this, we pick a single fiducial epoch in the color 
evolution and assess the performance of the color as a reddening estimator. Our polynomial 
fits to the color curves are very convenient for this purpose, since they allow us to interpolate 
reliable colors on a day-to-day basis over a wide range of epochs and explore which epoch 
gives the best results. 

If all SNe share the same intrinsic temperature at some epoch, we expect the subset 
of dereddened SNe to have nearly identical colors (Co); the remaining objects should show 
color excesses, E(C) = C — C , in direct proportion to their extinctions. A useful diagnostic 
to check our underlying assumption is the color-color plot: unreddened SNe should occupy 
a small region in this plane. If we further assume the same extinction law in the SN host 
galaxies, the subset of extinguished objects should describe a straight line originating from 
such a region. The figures of merit in this test are (1) the color dispersion displayed by the 
unreddened SNe, (2) the slope described by the reddened SNe (which is determined by the 
extinction law), and (3) the dispersion relative to the straight line (the smaller the better). 

We have identified four objects in our sample (SN 2003B, SN 2003bl, SN 2003bn, 
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SN 2003cn) consistent with zero reddening. These objects were selected for having (1) no 
significant Na I D interstellar lines in their spectra at the redshifts of their host galaxies, and 
(2) early-time spectra that show no signs of host-galaxy extinction. To judge the latter we 
used extinction val ues determined from fits of S N IIP atmosphere m odels to our early-time 



spectra as done by lDessart fc Hillierl ( 120061 ) and lDessart et al.l (120081 ) . Such models use the 
SN spectral lines to constrain the photospheric temperature and the continuum to restrict 
the amount of extinction. Examples of fits are shown in Figure HI 

We investigated two color-color plots (V — I versus V — R, and V — I versus B — V) 
over a wide range of epochs (from day —50 to —15) afte r corre cting the photometry for 



Galactic extinction and K -terms as described by lOlivared (120081 ). The best results were 
found for the V — I versus V — R diagram constructed from day —30 and shown in Figure [5J 
Fortunately, all SNe in our sample have enough data to interpolate color curves at day 30. 
At this epoch, approximately the end of the plateau, we obtain the linear behavior expected 
for a sample with the same intrinsic color but different degrees of extinction. Shown with 
a triangle is the single SN consistent with zero extinction which is, remarkably, one of the 
bluest objects in this diagram; the other three unextinguished SNe do not have i?-band 
photometry. A least-squares fit to the data yields a slope of 1.59 ± 0.09, which is lower than 
the E(V — I)/E(V — R) = 2.12 ratio expected for a Galactic extinction law (i?y = 3.1) 
shown as a vector in Figure |5j This suggests a somewhat different extinction law in the 
SN hosts compared to the Galaxy. The dispersion of 0.059 mag in V — I is a promising 
result as it translates into an uncertainty of A host (V) = 0.15 mag, which corresponds to the 
limiting precision of this method. The reduced x 2 °f 1-55 implies that the dispersion can 
be accounted almost solely by our error bars and that any instrinsic color dispersion in our 
sample is < 0.06 mag. The bottom line is that both the V — R and V — I colors fulfill the 
minimum requirements as reddening indicators. 

The best results from the B — V versus V — I analysis, obtained from day —30, are 
shown in Figure [6j The four SNe consistent with zero extinction, plotted with triangles, 
have weighted average colors of (B — V)q = 1.147 ± 0.053 mag and (V — J)o = 0.656 ± 0.053 
mag, where the dispersion was used as measure of the uncertainty. Note that there are five 
SNe in this diagram which are slightly bluer than the unreddened sample. The V — I color 
dispersion of 0.076 mag is greater than that obtained in Figure and is most likely due to 
the B — V color, since the B band is more sensitive to the metallicity of the SN, owing to 
many absorption lines that lie in this spectral region. Therefore, we believe that the greater 
dispersion in this diagram could be due to the different metallicities of our SN sample. A 
least-squares fit to the data yields a slope of 0.77 ± 0.04. This slope is quite different from 
the E(V — I)/E(B — V) = 1.38 ratio expected for the Galactic extinction law (shown as a 
vector in Figure [6]), in agreement with the departure from a standard reddening law seen in 
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the V — I versus V — R diagram. 

We conclude from our exploration that, while the B — V color is problematic, both the 
V — R and V — I colors offer a promising route for dereddening purposes. In what follows, 
we will employ solely the V — I color since only a small subset of our objects possess .R-band 
photometry. Although the evidence points to a non-Galactic reddening law, for now we will 
assume a standard reddening law (later we will relax this assumption; see § 14.3.21) . 

Using our library of SN II spectra, we computed synthetically the appropriate conversion 
factor between E(V—I) and Ay for SNe II and a standard reddening law [Ry = 3.1), yielding 



* = W^T) = 2 - 518 - (6) 

The host-galaxy extinction can be computed from A V (V — J) = f} v [{V — I) — (V — J) ] and, 
assuming an intrinsic color of (V — I)o = 0.656 ± 0.053 mag, its numerical expression with 
the corresponding uncertainty is 



A V (V - I) = 2.518 [(V-/) - 0.656], (7) 
a(A v ) = 2.518 J(?(y-i) + 0.053 2 + 0.059 2 , 



where V — I corresponds to the color of a given SN at day —30 (corrected for i^-terms 
and foreground extinction) and <J(v~i) combines the in strumental err ors in the V and I 



magnitudes, the uncertainty in Aq and the ii'-terms (see IOlivaresll2008l . Chapters 3.1.1 and 
3.1.2, respectively), and the uncertainty in tpx (§ I3.1.ip . The uncertainty in the intrinsic 
(V — I)o color was also included in the error of Ay along with the V — I standard deviation 
in the V — I versus V — R diagram. 

The host-galaxy reddening values obtained from this technique are listed in Column 5 
of Table H] (with the uncertainties given in parentheses for the whole sample of 37 SNe). 
There are eight SNe with V — I colors bluer than (V — J)q which stand out in this table 
for their inferred negative reddenings. Although not physically meaningful, these negative 
values are statistically consistent with zero or moderate reddenings. In fact, seven of these 
eight objects differ by < 1.1a from zero reddening. Only SN 1992af differs by 1.8<r from 
A host = mag. 



As a comparison, iKrisciunas et al.l (120091 ) employed a similar procedure to determine 
the extinction to SN 20031m. While we get Ay = 0.46±0.21, they obtained Ay = 0.54±0.14, 
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values that are consistent within the uncertainties. In the following section we compare this 
method with other dereddening techniques. 



4. ANALYSIS 
4.1. Comparing Dereddening Techniques 



SN IIP atmosphere models w ere fitted to our library of spectra as done by lDessart fc Hillier 



(120061 ) and iDessart et all (120081 ). In these fits the spectral lines are used to constrain the 
photospheric temperature and the corresponding continuum is employed to estimate the ex- 
tinction. Thus, the quality of the fit depends on (1) the determination of the photospheric 
temperature by means of the spectral lines, and (2) how well the continuum is represented 
by the extinguished blackbody emission. A crucial condition for this technique to work is 
the spectrophotometric quality of the spectra. 

I n general, our obs ervations were obtained with the slit oriented along the parallactic 



angle (jFilippenkol Il982h and the relative shape of the spectra should be accurate. However, 
this was not always possible; moreover, sometimes the spectra were contaminated by light 
from the host galaxy. Accordingly, we first checked the flux calibration of each spectrum 
by synthesizing magnitudes and comparing them to the observed magnitudes, duly interpo- 
lated to the time of the spectroscopy. We typically found an agreement better than 0.05 mag 
between the observed and synthetic colors, thus confirming our confidence in the flux cali- 
bration. In the 64 cases (18% of all spectra) where we found significant differences between 
synthetic and observed colors (> 0.1 mag for any color), we applied a low-order polynomial 
correction to the spectrum. Basically, this "mangling" correction used the observed photom- 
etry to change the slope of a spectrum. After checking the flux calibration (and mangling it, 
if needed), the next step was to correct for Galactic absorption and deredshift the spectra. 
This database was then used for the atmosphere model fits. Examples of the fits are shown 
in Figure HI 

The result ing spectroscopic redden ings are summarized in Column 2 of Table HI As 



pointed out by IDessart fc Hillier! ( 120061 ). the spectrum-fitting technique works much better 
with early-time spectra than with late-time observations. At late times the photosphere 
has receded in mass, exposing inner and more metal-rich layers, which translates into an 
overabundance of heavy-element absorption lines. Therefore, the fitting to the continuum 
(essentially a temperature fitting) is hampered by the increased line opacity at late times. We 
thus divided our sample into four subcategories based on the epoch and the flux quality of the 
spectra used in the reddening determination: (a) gold: early-time spectra, without mangling 
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correction; (b) silver: early-time spectra, with mangling correction; (c) bronze: late-time 
spectra, without mangling correction; and (d) coal: late-time spectra, with mangling cor- 
rection. Note that the main criterion is whether the spectrum was obtained at early or late 
times, and the second criterion corresponds to the flux-calibration quality. We consider the 
unmangled spectra to be of higher quality than the mangled ones, because the former do 
not require any corrections and suggest that the observations were better performed. 

The above sub classification is given for each SN in Co lumn 3 of Tab l e HI Note that 



the uncertainty is directly related to this subclassification. iDessart et al.l (120081 ) estimate 
an error of o~e(b-v) = 0.05 mag when using early-time spectra, and (Je(b-v) = 0.10 mag 
when using late-time spectra. We refine this argument, ramping up from 0.05 to 0.10 mag, 
depending on the number of spectra employed for each extinction determination. These 
values are given in Table S] in units of Ay = 3.1 E(B — V). 

Figure [7] shows a comparison between the spectroscopic reddenings Ay (spec) and our 
color reddenings Ay(V — I) derived in § I3.2[ for the 17 SNe belonging to the top three 
subclasses (gold + silver + bronze). Good agreement is found between both techniques, with 
a dispersion of 0.38 mag. The resulting xl — 1-4 suggests that this dispersion is consistent 
with the combined errors of both techniques. The exceptions are two objects, SN 1999br and 
SN 2003ho. SN 1999br has only plateau-phase photometry, making it hard to determine tpx- 
Although the uncertainty in tpp is quite large (20 days), this does not have a great impact 
on the error of the V — I color due to the flatness of the color curve at these epochs (+0.0024 
mag per day). Another possible reason for the disagreement could be an intrinsically redder 
V — I color than that of the bulk of the SNe, caused by the extreme properties (low luminosity 
and expansion velocity) of SN 1999br. The second discrepant object (SN 2003ho) belongs to 
the bronze group, so it is possible that the difference could be due to the use of a late-time 
spectrum in the determination of the spectroscopic reddening. 

Another way to estimate host-galaxy reddening is from the Na I D AA5889, 5896 inter- 
stellar absorption doublet observed in the SN spectrum at the host-galaxy redshift. When- 
ever the line was detected we measured its equivalent width; in those cases where we did not 
detect the Na I D line we assigned it a null value. In all cases we estimated the uncertainty 
in the equivalent width (EW) based on the S/N of the continuum around this line. We 
converted the s e mea surements into visual extinctions, Ay(N& I D), using the calibration of 



Barbon et al.l fll990h : 



E(B-V) « 0.25 EW(Na I D), 



where the units of E(B — V) and EW are mag and A, respectively. Although iTuratto et al 
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( 120031 ) showed that this relation is bimodal, and lBlondin et al.l (120091 ) found that the E(B — 
V) estimates based on EW(Na I D) were highly inaccurate, the E(B — V) values derived 
from Equation [S] (tabulated in Column 4 of Table HJ) provide a rough estimate of E(B — V) 
to compare with our results. The comparison between Ay(N& I D) and our technique, shown 
in Figure [HJ exhibits a dispersion of 0.53 mag (Xv = 3.8), which is much higher than the 
a = 0.38 mag obtained from the previous comparison. The disagreement can be attributed 
to the fact that the absorp tion line traces the gas c ontent along the line of sight, but does 
not necessarily probe dust (IMunari fe Zwitterlll9971 ). Furthermore, reddenings derived from 
the EW of the Na I D lines measured from our low-dispersion spectra simply cannot be 
expected to b e accurate; the D lines produced by a typical interstellar cloud are saturated 
(IHobbs Ill974l ). The only way that one can hope to derive a reliable reddening from the D 
lines is via very high-dispersion spectra that resolve the lines and allow the column density 
to be determined. 



4.2. Luminosity versus Expansion Velocity 



After developing the dereddening met hod based on late - time colors, we can now revisit 
the LEV relation originally discovered by iHamuy fe Pintol ( 120021 ). which is at the core of 



the SCM. For this purpose we applied AKA corrections to our photometry (jOlivaresI 12008 
Chapter 3.1), we used our analytic fits (§ 13. 1.1)) to interpolate BVI magnitudes on day —30, 
and we employed the cosmic microwave background (CMB) redshifta^l in Table [5] to convert 
apparent magnitudes to absolute values (assuming Hq = 70 km s -1 Mpc -1 from our final 
result). In order to assess the dereddening technique, we preserved the relative reddening 
estimations by keeping the negative values of host-galaxy absorption. We performed a power- 
law fit to interpolate a velocity contemporaneous with the photometry (day —30) as described 
in § 13.1.11 Only for SN 1999cr did we have to use an extrapolated velocity (see top-right 
panel in Figure [3]). From our original sample of 37 SNe IIP we were able to use 30 SNe to 
build this relation, since five of them do not have measured Fe II velocities at day —30, and 
two others have extremely low redshifts (czcmb < 300 km s _1 ). 

Figure [9] shows the resulting LEV relation (absolute m agnitude versus expansi on veloc- 
ity) for all BVI bands. Evidently we recover the result of IHamuy fe Pintol (120021 ). namely 
that the most luminous SNe have greater expansion velocities. In their case the data were 



10 The CMB redshifts were computed by adding the heliocentric redshifts given in Tableland the projection 
of the velocity of the Sun relative to the CMB in the direction of the SN host galaxy. For the latt er we 
adopted a velocity of 371 km s _1 in the direction (l,b) — (264.1°, 48.3°) given by Fixsen et al. ( 19961) plus 
the uncertainty of the Local Group velocity (187 km s _1 ). 
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modeled with a linear function. Our sample suggests that the relation may be quadratic, but 
we need more SNe at low expansion velocities to confirm this suspicion. Linear least-squares 
fits to our BVI data yield the following solutions: 



M(£)_3o-51og(# /70) 
M(V)_30-51og(F /70) 
M(/)_ 30 -51og(iV70) 



3.50(±0.30) log t; Fcn (-30)/5000 - 16.01(±0.20), (9) 
3.08(±0.25) log [v Fe ii(-30)/5000] - 17.06(±0.14), (10) 
2.62(±0.21) log [-u Fen (-30)/5000l - 17.61(±0.10), (11) 



which are shown with solid lines in Figure |9j The relation found by iHamuy fc Pintol (120021 ) 
for the V band is indic ated with the da s hed li ne in the middle panel of the same figure. 
Given that the study of IHamuy fc Pintol ( 120021 ) was performed using data at day 50 after 
explosion (around day —60 on our time scale), it is not unexpected that their LEV relation 
is shifted to higher expansion velocities. Some of the difference in slope is explained by the 
inclusion of SN 2003bl in our sample, which flattens the correlation. This relation exhibits a 
dispersion of 0.3 mag, similar to that reported before. Such low scatter is very encouraging 
as it implies that the expansion velocities can be used to predict the SN luminosities, to 
standardize them, and to derive distances. 



Recently, iKasen fc Woosleyl ( 120091 ) have comp uted hydrodynam i cal m odels for SNe IIP, 
which support the LEV relation initially reported by lHamuy fc Pintol (120021 ) . Also, M. Bersten 
(2010, in preparation) has used her own hydrodynamical models to investigate the LEV re- 
lationship. By varying the explosion energy, she has obtained an LEV relation with a slope 
consistent with that reported here. 



4.3. Hubble Diagrams 

In order to assess the performance of the SCM, we employ (a) host-galaxy redshifts in 
the CMB frame with an uncertainty of 187 km s _1 , (b) BVI apparent magnitudes and V — I 
colors corrected for K-terms and Galactic reddening, and (c) Fe II velocities (in units of 
km s _1 ). Table [5] lists these values for the 37 SNe of our sample, of which 35 meet an initial 
requirement of czqmb > 300 km s _1 . The weighted dispersion of the data points around a 
linear fit in the HDs allows us to evaluate the precision of different approaches. 
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4-3.1. Using A V (V - I) and A v (spec) 

The top-left panel of Figure [10] shows a HD constructed from 5-band magnitudes in- 
terpolated to day —30. The top-right panel shows magnitudes additionally corrected for the 
LEV relation, and the bottom-left panel includes further corrections for host-galaxy extinc- 
tion (using the V — I color calibration given in § I3.2p . In each case we perform a weighted 
linear least-squares fit of the form 



m + a log(v FeII /5000) - A host = 5\ogcz + zp, (12) 

where m is the apparent magnitude corrected for K-terms and Galactic absorption, v-p e n is 
the expansion velocity, A host is the host-galaxy absorption, and z is the CMB redshift. Note 
that the SN redshift enters logarithmically in this model, implying that low-redshift SNe 
carry less weight in the fit. The only fitting parameters are a and the zero point (zp), but 
in the top-left panel we set a = and we fit only for the zp. The significant decrease in 
the dispersion, from 0.72 to 0.36 mag, clearly demonstrates the need to add the velocity and 
host-galaxy extinction terms. 

An inspection of the V-band diagram (Figure fTTj) shows a large scatter of 0.54 mag in 
the top-left panel. When we correct for expansion velocities, the scatter drops to 0.50 mag. 
This is certainly not unexpected given the LEV relation reported in the previous section. It 
is encouraging to notice that the dispersion drops from 0.50 to 0.45 mag when we include our 
host-galaxy extinction corrections; this indicates the utility of our color-based dereddening 
technique. The reduced \ 2 value of 2.45 implies that most of the scatter is caused by the 
observational errors. We performed the same analysis using other epochs and found that 
day —30 yielded the lowest dispersion. 

If we turn our attention to the I band (Figure [T2l) . the final dispersion is 0.45 mag, iden- 
tical to that found in the V band. It seems that the dispersion could have some dependence 
on wavelength, since it decreases from 0.45 in VI to 0.36 mag in B. The scatter of ~ 0.4 mag 
in BVI is comparable to, but somewhat larger than, the ~ 0.34 m ag dispersion obtained in 



previous SCM studies (IHamuy fc Pintdl2002l ; lHamuyil2003bl . 120051 ). It is important to notice 



that the dispersions are independent of the intrinsic (V — I)q color calculated in § 13.21 

In the bottom-right panels of Figures [TQl [HI and [Tj2] we examine the SCM using the 
spectroscopic reddenings determined for a set of 28 SNe IIP from Table H] instead of our color- 
based extinctions. The resulting dispersions in BVI are (0.67, 0.54, 0.42 mag), compared 
with (0.36, 0.45, 0.45 mag) when using color-based extinctions. All of the fitting parameters 
derived from both dereddening techniques are compiled in Table EJ 
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4-3.2. Leaving Ry as a Free Parameter 

The previous analysis suggests that the dispersions are somewhat larger than the ob- 
servational errors. One possible source of scatter could be the reddening law. As shown in 
§ !3.2l and Figure 0, we have evidence for a somewhat different extinction law in the SN hosts 
compared to the Galaxy. Here we take this idea a step further and analyze the HD leaving 
Ry as a free parameter. To accomplish this we model the data with the expression 



m + a log (^Feii/5000) — {3(V — I) = 5\ogcz + zp. 



(13) 



Note that we are replacing v4h os t in Equation ffl2|) with the term (3(V — I). The /3 factor is 
a new free parameter to be marginalized along with a and zp, and V — / is the color on 
day — 3 unc orrected for host-galaxy dust absorption. This is the same approach adopted by 
Trippl ( 119981 ) for SNe la but using decline rate instead of expansion photospheric velocity. 
The underlying assumption in our case is that the absolute magnitudes of SNe IIP are a 
two-family function of Fe II expansion velocity and V — I color. The latter is interpreted as 
a color excess which varies according to dust in the host galaxy. Once (3 is known we can 
use it to solve for host-galaxy extinction from A(X) = j3 (V — I) — 0.656 
the intrinsic V — I color of SNe IIP (see Equations (jSJ) and ([7]) in § 13.21) . 



where 0.656 is 



Figure [13] shows the BVI HDs for our set of 30 SNe IIP. We obtain dispersions of 
(0.28, 0.31, 0.32 mag) in BVI, respectively, compared with (0.36, 0.45, 0.45 mag) when R v 
is kept fixed. The xl values increase from (1.6, 2.5, 2.5) to (2.8, 5.4, 8.7). This change is 
due to the fact that we are not using the intrinsic color, which gives the major contribution 
to the errors, even more in the I band. Although we expect a reduction in the scatter due 
to the inclusion of an additional parameter, the large drop in the dispersion is remarkable. 
The last three lines of Table [6] show the parameters we obtain by minimizing the dispersion 
in the HD for each band. When we restrict the sample to objects with czqmb > 3000 km s _1 
(leaving aside SNe with potentially greater peculiar velocities), we end up with 19 SNe in 
the B band and 20 SNe in the VI bands. The resulting HDs in BVI show dispersions of 
(0.25, 0.28, 0.30 mag), respectively, and the zero points of the HDs change by only 0.02 mag. 
Not surprisingly, these are lower than the (0.28, 0.31, 0.32 mag) dispersions obtained from 
the whole sample. 

Efforts to remove the peculiar velocities of low-redshi f t gala xies were made by means 
of the parametric model for peculiar flows of iTonry et al.l ( 120001 ). The ingredients of this 
model are (1) the observed recession velocity of the galaxy, (2) t he Hu bble constant, and 
(3) the surface-brightness- fluctuations distance fro m ITonry et al.l (120011 ) . The output is an 



estimate of the peculiar velocity of the galaxy (see iHamuyl l2003al . § 1.4 for details). After 
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correcting 10 host galaxies in our sample, the results are not conclusive. The scatter in the 
HDs decreased only 0.01 mag for the / band. 

Our assumption is that the term /3(V — I) in Equation f[T3~j) corresponds to the extinction 
in a broad-band magnitude (with central wavelength A). Thus, /3 is the ratio A(X)/ E(V — I) 
(see Equation ([6]) in § 13.21) and is related to the shape of the extinction law. For each of the 
BVI bands, we used our library of SN II spectra to compute synthetically the value of (3\ as 
a function of Ry = A V /E(B — V); see Figure dH This allowed us to convert the (3\ values 
resulting from our fits into R v values. Our fits yield (3 X = (2.67±0.13, 1.67±0.10, 0.60±0.09) 
for the BVI bands, respectively, which translate into R v = (1.3Stoll,lAitoslA-^-oil)- 
These values are remarkably consistent with each other and significantly lower than the 
Rv = 3.1 value of the standard reddening law. Independent evidence for a low- Ry law has 
already been reported from studies of SNe la (see § E]). 

The dispersion of 0.28-0.32 mag in the HDs translates into a precision of 13-14% in the 
determination of relative distances. However, as noticed after inspection of Figure [151 the 
magnitude residuals are highly correlated among BVI and generally decrease with increasing 
redshift. The likely explanation is that we are seeing the effects of peculiar velocities of the 
host galaxies, in which case the precision of the SCM could be better than the 13-14% 
derived from the scatter in the HDs. To examine this point, for each SN we took the czqmb 
values from Table and calculated the distance modulus residual 



A/z = m + a log (-UFeii/5000) — f3(V — I) — 5 log cz — zp, 



(14) 



for all BVI bands. Figure [16] shows strong correlations of these magnitude differences 
against each other. The root-mean square (rms) scatter in the Afi(B) vs. Afi(V) relation 
is 0.12 mag, and for Afi(B) vs. A/i(J) it is 0.19 mag. This translates to 6% and 9% in 
distance, respectively. If this correlation of the residuals is due to peculiar velocities, then 
the true distance precision of SCM is 6-9%, and not 13-14% as derived from the Hubble 
diagrams (which include scatter produced by the peculiar velocities). Th i s result is very close 
to the ~ 7- 10% precision of methods using SNe la (e.g., |Phillipslll993t lHamuy et al.l Il996 



Riess et al.l Il996l ; iPhillips et al.l Il999[ ) when they were first used to reveal the acceleration 
of the Universe, which demonstrates that SNe IIP have great potential for determining 
extragalactic distances and therefore cosmological parameters. 
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4.4. The Hubble Constant 



The Hubble constant is a parameter of central importance in cosmology which can be 
determined from our HDs. This can be accomplished as long as we can convert apparent 
magnitudes into distances. The calibration was done with two objec ts for which we found 
Cepheid distances in the literature: SN 1999em {fx = 30.34 ±0 .19 mag; iLeonard et al.ll2002af ) 



and SN 2004dj {fx = 27.48 ± 0.24 mag; Freedman et al J BOO If ). For each calibrating SN we 
can solve for H using the expression 



log Hq 



0.2 



m + a log (t; Fen /5000) - 0{V - I) - fx - zp + 25 



-30 



(15) 



where m is the apparent magnitude, -upeii is the expansion velocity, and V — I is the color 
of the calibrating SN, all measured on day —30 and given in Table [5J a, (3, and zp are the 
fitting parameters given in the bottom three lines of Table [6j 

Table [7J summarizes our calculations for the two calibrators. We see that the corrected 
BVI absolute magnitudes of SN 1999em and SN 2004dj differ from each other by 0.86- 
1.13 mag. This is considerably greater than the scatter of 0.3 mag in the LEV relation, 
but statistically not impossible; see Figure [151 where we plot with triangles the corrected 
absolute magnitudes of the two calibrating SNe, on top of the whole sample of SNe employed 
in the HDs. This difference leads to H Q values in the range 62-105 km s _1 Mpc -1 . With 
only two calibrating SNe, the SCM still has much room for delivering a more accurate and 
precise value for H . 



4.5. Distances 



We are in a position to calculate distances to all of the SNe of our sample. This can be 
done with the expression 



fx = m + a log (i; Fen /5000) - [3{V -I)- (M corr ), (16) 

where (M corr ) is the corrected absolute magnitude of SNe IIP on day —30. For this purpose 
we employ the weighted mean of the two calibrating SNe given in Table [71 The resulting 
values are given in Tabled for 32 SNe IIP; the last column shows the weighted average of 
their distance moduli and its uncertainty, assuming no covariance. As Equation [16] shows, 
it is straightforward to shift the resulting distance moduli to accommodate a new value for 



-22 - 



(M corr ). Unfortunately, we are not able to estimate expansion velocities for five objects, and 
consequently not their distances as well. 

As mentioned in § [TJ we can evaluate the accuracy of the SCM from a comparison 
with EPM distan ces. For this purpose we employ EPM distances recently calculated by 
Jones et al.l (120091 ). which are summarize d in Table [9] along w ith our distances for the 11 ob- 
jects in common between SCM and EPM. I Jones et al.l (120091 ) determine EPM distances with 
two different sets of atmosphere models, both developed based on an analytic approach to 
calculate dilution factors as a function of a color temperature in order to correct the assumed 
blackbody emission. Their EPM distances were derived from average dilution factors and 
n ot from specific models for each supernova. The EPM distances obtained using the models 
of Eastman et al.l (119961 ) (E96; Column 2) are 12 ± 5% shorter than our SCM distances. On 
the other hand, the EPM distances determined using the models of iDessart fe Hillierl (120051 ) 
(D05; Column 3) are 40 ± 10% greater than our SCM distances. Both percentage shifts are 
calculated weighted by the uncertainty in the EPM and SCM distances. These systematic 
differences can be seen in the top panel of Figure [T7J and more clearly in the fractional 
differences (g?epm — ^scm)/^scm plotted in the bottom panel. 

The systematic differences among the two sets of EPM distances can be solely attributed 
to the atmosphere models of E96 and D05. New radiative transfer models of SNe IIP are 
necessary to understand the origin of this discrepancy. Besides the systematic errors in 
both EPM implementations, we can gain an understanding of the internal precision of EPM 
and SCM after removing the systematic differences and bringing all distances to a common 
scale. For this purpose we correct the EPM distances to the SCM distance scale by removing 
the percentage shifts. The comparison is shown in Figure [181 t ne distance differences have 
dispersions of ~ 13% and ~ 16% using D05 and E96, respectively. This implies that SCM and 
EPM produce relative distances with an internal precision < 13-16%. This agrees with the 
dispersions of 0.25-0.30 mag seen in the HDs restricted to SNe with c^cmb > 3000 km s _1 . 

Further comparison betwee n EPM and SCM can be made by taking the distance deter- 
mination to SN 1999em. While IDessart fc Hillierl (120061 ) give a value of 11.5 ± 1.0 Mpc, our 
SCM returns 12.6 ± 0.6 Mpc. This agreement has a significance better than la, taking into 
account both sources of error. 
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DISCUSSION 



5.1. Variations of the Extinction Law 



1.4±0.1) turns 



The reddening law that minimizes the dispersion in the BVI HDs (Ry 
out t o be very di fferent from the s tanda rd Galactic extinction law (Ry = 3.1; iCardelli et al. 



19891 ). Recently, iPoznanski et al.l ( 120091 ) carried out a similar analysi s using SNe IIP, which 
confirmed the low value obtained for Ry presented bv lOlivaresi (|2008h. There hay e also been 
several studies on this subject using SNe la. Most recently, iFolatelli et al.l ( I2010l ) solved for 
Ry in the same manner as our approach (i.e., by minimizing the dispersion in the HD) and 
obtained a value of Ry ~ 1.5, remarkably consistent with our value. Further e vidence for low 



R\ values using SNe la were reported by others applying a similar procedure ( Phillips et al. 
19991 ; lAltavilla etaDbood ; iReindl et allbood : Ifficken et allhood : IWang et al.ll2009h . 



Other studies of dust reddening based on diverse methods were performed for nearby 
galax i es giving values for Rb ranging between 2.4 and 4.3 faifattJl990 l: iDella Valle & Pana gia 
19921 ; iBrosch fc Loingerl Il99ll ; iFinkelman et al.l 120081 ) . The ratio of total to selective ab- 
sorption varies significant ly over the range Ry ~ 2.6 — 5.5 even within our own Galaxy 
jciayton fc Cardellilll988h . 



A lower value of Ry could be due to dust grains smaller than those in our Galaxy, since 
for a given va lue of Ay the E(B — V) reddening increases if the size of the grains decreases 
( )Drainell2003l ). However, it would be very unlikely that our Galaxy is so special. 



Most probably we are dealing with normal dust grains in the circumstellar medium, and 
efore multiple scattering of photons makes the reddening law steeper and Ry smaller 



(Wang 


2005; 


Goobar 


2008) 



be further tested with photometry at such wavelengths. 



5.2. Internal Precision of the SCM 

The BVI HDs have a scatter of ~ 0.3 mag, implying a precision in individual distances 
of ~ 15%. Part of this scatter might be due to the peculiar velocities of the SN host 
galaxies, and the intrinsic precision of SCM could be even better. In fact, when we perform 
a comparison for SNe in common between SCM and EPM, the distance differences are 13- 
16% (after removing the systematic difference among the SCM and EPM). This comparison 
is independent of the SN host-galaxy redshifts and implies that the internal precision in 
either of these two techniques must be lower than 13-16%, since these differences comprise 
the combined uncertainties of both methods. This is an encouraging result; it implies that 



-24- 



both EPM and SCM can produce high-precision relative distances, thus offering a new route 
to cosmological parameters. 

5.3. The LEV Relation 

We note that the dispersion in the HDs increases with wavelength. This seems con- 
trary to expectations given that (1) the extinction effects decrease with wavelength, and (2) 
metallicity affects the B band more than the other optical bands. Perhaps the luminosity is 
a function of not only velocity but also metallicity: My = f(v, [Fe/H]). If so, the velocity 
term in Equations ( Tl2|) and ( [TBI might be removing metallicity effects more efficiently at 
shorter wavelengths. This could be tested with [Fe/H] measurements of SNe IIP. On the 
other hand, hydrodynamical simulations of SNe IIP suggest that there is an intrinsic scatter 
in the relation, which does not seem to be driven by any parameter of the model. Given 
the evidence we have, we can only claim a dependence of I on t;. The obvious physical 
explanation is that the internal energy is correlated with the kinetic energy, implying that 
the ratio of the internal energy to the kinetic energy is approximately independent of the 
explosion energy. As SN 1999br may suggest, underluminous events might (1) change the 
shape of the LEV relation, or (2) require a different explosion mechanism. 

6. CONCLUSIONS 

We established a library of 196 SN II optical spectra and developed a code which allows 
us to correct the observed photometry for Galactic extinction, i^-terms, and host-galaxy 
extinction. We developed fitting procedures to the light curves, color curves, and velocity 
curves which allow us to precisely determine the transition time between the plateau and 
the tail phases. The use of this parameter as the time origin permitted us to line up the SNe 
to a common phase. The additional benefit of these fits is the interpolation of magnitudes, 
colors, and velocities over a wide range of epochs. The methodology explained above yields 
the following conclusions. 

1. The comparison between our color-based dereddening technique and the spectroscopic 
reddenings is satisfactory within the uncertainties of both techniques. This is particu- 
larly encouraging since our method uses late-time photometric information, while the 
other method uses completely independent early-time spectroscopic data. 



2. Using our new sample of SNe we rec over the luminosi t y vers us expansion- velocity trend 
(LEV relation) previously found by lHamuy &: Pintol (120021 ) . 



-25 - 



We construct BVI Hubble diagrams using two sets of host-galaxy reddenings. We 
demonstrate that the extinction corrections based on V — / color do a better job than 
the spectroscopic determinations, reaching dispersions of ~ 0.4 mag in the Hubble 
diagrams. This scat ter is somewhat higher than the 0.35 mag found previously by 



Hamuyl (j2003bl . 120051 ). A much smaller dispersion of 0.3 mag is achieved when we use 
V — I colors to estimate reddening and allow Ry to vary. We obtain Ry = 1.4 ± 0.1, 
much smaller than the Galactic value of 3.1. The low value of Ry can be explained by 
a different nature of the dust grain s in the host galaxies of SN e IIP or by scattering of 
light by circumstellar dust clouds (I Wang 1120051 ; lGoobarll2008l ). 



A clear correlation is found among residuals of BVI HDs. After removing the corre- 
lation, the rms in the HDs decreases to 6-9% in distance, comparable to the precision 
achieved with SNe la. 

After calibrating our Hubble diagrams with Cepheid distances to SN 1999em and 

SN 2004dj, our BVI photometry leads to a Hubble constant in the range 62-105 km s _1 Mpc -1 ; 

this agrees with most of the values obtained by means of diverse modern techniques. 

Finally, we calculate the distance modu li to our SNe, and make a comparison against 
EPM distances from I Jones et al.l ( 120091 ). The 11 SNe in common show a systematic 
difference in distance between EPM and SCM, depending on the atmosphere models 
employed by EPM. Correcting for these shifts we bring the EPM distances to the SCM 
distance scale, from which we measure a dispersion of 13-16%. This spread reflects 
the combined internal precision of EPM and SCM. Therefore, the internal precision in 
either of these two techniques must be < 13-16%. 



This analysis confirms the utility of SNe IIP as cosmological probes, providing strong 
encouragement to future high-redshift studies. We find that one can determine relative 
distances from SNe IIP with a precision of 15% or better. This uncertainty could be further 
reduced by including more SNe in the Hubble flow. In its current form the SCM requires 
both photometric and spectroscopic data. Since the latter are expensive to obtain (especially 
at high redshifts), it would be desirable to look for a photometric observable as a luminosity 
indicator instead of the expansion velocities. 



We thank Luc Dessart, who generously provided the reddening estimates calculated by 
means of spectral fitting and the models plotted in Figure HI F.O.E. also thanks him for 
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Fig. 1. — BVRI light curves of SN 1999em. BVRI magnitudes are respectively shown in 
blue (top- left), green (top- right), brown (bottom- left), and red (bottom-right). The colored 
solid line corresponds to the analytic fit. The black lines represent the decomposition of the 
main function into its three components: the Fermi-Dirac function (dotted line), a straight 
line (long-dashed line), and a Gaussian function (short-dashed line). The origin of the time 
axis corresponds to tp? in Julian days for SN 1999em. Since the instrumental error bars are 
smaller than the symbols, they are not shown. 



- 33 - 




I i i i I i i i I i i i I i i i L 

-80 -60 -40 -20 

JD-Up [days] 



Fig. 2.— B- V, V-R, and V-I color curves of SN 1999em (filled circles), SN 1991al (open 
circles), and SN 2003hl (filled triangles) corrected for the Aq and i^-terms. This comparison 
demonstrates that each SN displays a different color evolution, precluding the determination 
of color excesses from a simple color offset after correcting for the Aq and K-terms. For 
SN 1991al the error bars are always smaller than the symbols and therefore not plotted. 
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Fig. 3. — Four expansion velocity curves representative of the SN sample measured by means 
of the Fe II A5169 line profile. While solid lines correspond to a power-law fit (equation [5]) , 
dashed lines show the transition time. The to parameter of the fit was choosen as the time 
origin due to plotting restrictions. For SN 1992af (top-left panel) we fit only two velocities 
(see § 13.1.31 for more details). The top-right panel shows a 15-day extrapolation beyond 
the last data point for SN 1999cr. The bottom-left panel shows a 10-day extrapolation for 
SN 2003gd prior to the first data point. SN 1999em velocities are shown in the bottom-right 
panel as an example of good sampling. 
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Fig. 4. — Type IIP SN atmosphere models (wide red line) fitted to our spectra (black line) 
in the rest frame. The top panel shows an example of an unsatisfactory fit to one of our late- 
time spectra. The middle panel shows a better fit to an early-time spectrum of SN 1991al. A 
much better fit is achieved for an early-time spectrum of SN 2003bl as shown in the bottom 
panel. Strong telluric lines are marked. 
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Fig. 5. — V — I versus V — R diagram for 21 SNe IIP having VRI photometry corrected 
for the Aq and i^-terms. The blue line is a least-squares fit to the data, with a slope of 
1.59 ± 0.09. The blue arrow has a slope of 2.12 and corresponds to the reddening vector for 
Ay = 1.0 mag and a standard Galactic extinction law {Ry = 3.1). The blue triangle shows 
the one SN of this subsample consistent with zero extinction having i?-band photometry. 
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Fig. 6. — V — I versus B — V diagram for 29 SNe IIP having BVI photometry corrected 
for the Aq and i^-terms. The red line is a least-squares fit to the data, with a slope of 
0.77 ± 0.04. The red arrow has a slope of 1.38 and corresponds to the reddening vector for 
Ay = 1.0 mag and a standard Galactic extinction law (Ry = 3.1). The red triangles show 
the four SNe consistent with zero host-galaxy extinction. 
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Fig. 7. — Comparison between two dereddening techniques via Ay (mag): the spectrum- 
fitting method (ordinate) and our technique (abscissa; § 13. 2p . The figure shows three sub- 
classes defined by the quality of the spectroscopic data used: the gold, silver, and bronze 
samples (see § 14. ip . The solid black line is fixed to show a slope of unity. 
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Fig. 8. — Comparison between two dereddening techniques via Ay (mag): the Na I D 
interstellar line (ordinate) and our technique (abscissa; § 13.21) . In each case, the equivalent 
width of the Na I D interstellar line is tr ansformed into a visual extinction according to the 
calibration given by iBarbon et al.l (ll990f ) . 
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Fig. 9. — BVI-b&nd absolute magnitudes (ordinate) against the Fe II expansion velocity 
(abscissa) for 29-30 SNe. We use magnitudes corrected for Galactic extinction, K-terms, 
and host-galaxy extinction together with a value of the Hubble constant of 70 km s -1 Mpc -1 . 
The dashed line in the middle panel represents the LEV relation for the V band found by 



Hamuy fe Pintol (120021 ). obtained from magnitudes and velocities measured at day 50 past 



the explosion (approximately day —60 on our own time scale). 
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Fig. 10. — B-band Hubble diagrams. The top- left panel shows magnitudes interpolated 
to day —30 previously corrected for Galactic extinction and X-terms; the top-right panel 
shows magnitudes additionally corrected for expansion velocities; the bottom-left panel in- 
cludes further corrections for color-based host-galaxy extinction (using the V — I color cal- 
ibration given in § 13. 2j) ; and the bottom-right panel replaces the color-based extinction for 
the spectroscopic reddenings. The horizontal bars in each panel correspond to an adopted 
uncertainty of 187 km s _1 in the CMB redshift. 
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Fig. 11. — V-band Hubble diagrams. The top- left panel shows magnitudes interpolated 
to day —30 previously corrected for Galactic extinction and X-terms; the top-right panel 
shows magnitudes additionally corrected for expansion velocities; the bottom-left panel in- 
cludes further corrections for color-based host-galaxy extinction (using the V — I color cal- 
ibration given in § 13. 2j) ; and the bottom-right panel replaces the color-based extinction for 
the spectroscopic reddenings. The horizontal bars in each panel correspond to an adopted 
uncertainty of 187 km s _1 in the CMB redshift. 
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Fig. 12. — /-band Hubble diagrams. The top-left panel shows magnitudes interpolated to 
day —30 previously corrected for Galactic extinction and i^-terms; the top-right panel shows 
magnitudes additionally corrected for expansion velocities; the bottom-left panel includes 
further corrections for color-based host-galaxy extinction (using the V — I color calibration 
given in § !3.2p ; and the bottom-right panel replaces the color-based extinction for the spectro- 
scopic reddenings. The horizontal bars in each panel correspond to an adopted uncertainty 
of 187 km s" 1 in the CMB redshift. 
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Fig. 13. — BVI Hubble diagrams leaving Ry as a free parameter. In the top, middle, 
and bottom panels we show the HDs using BVI-band magnitudes of 29, 30, and 30 SNe, 
respectively. The horizontal bars in each panel correspond to an adopted uncertainty of 
187 km s" 1 in the CMB redshift. 
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Fig. 14. — Ry versus (3. The solid lines are computed from our library of SN II spectra for 
each of the BVI bands. The /3 parameters derived from the Hubble diagrams in Figure [12] 
are shown with red open circles for each of the BVI bands. The red dotted line shows the 
Ry- weighted mean for the BVI values of (3. The values for the standard Galactic extinction 
law are also shown with green dots. 
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Fig. 15. — BV /-corrected absolute magnitudes (ordinate) against the Fe II expansion veloc- 
ity (abscissa) for 29-30 SNe (black dots). We use the H values of Table Hand the CMB 
redshifts to compute the distance moduli. The dashed lines show the mean corrected ab- 
solute magnitude for the black data points. The red triangles are the two calibrating SNe, 
whose corrected absolute magnitudes were calculated using their Cepheid distances (see Ta- 
ble [7]). Note that, as expected, the two calibrating SNe fall on each side of the corrected 
absolute magnitude distributions. 
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Fig. 16. — Residuals from the BVI HDs. The correlation between the V-filter and /-filter 
residuals against the S-filter residuals is shown in the upper and lower panels, respectively. 
The red lines represent the fit from which the scatter is measured. 
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Fig. 17. — Direct comparison between SCM and EPM distances in Mpc calculated by 



Jones et al.l ( 120091 ). In magenta are shown the EPM distances computed from the atmo- 
sphere models of E96, while the green triangles are EPM distances obtained from the D05 
models. The bottom panel shows the fractional differences between both techniques against 
log(dscM)- In both panels the dotted lines trace the systematic shifts of the EPM distances 
with respect to the SCM distances, ~ 40% and ~ 12% using D05 and E96 atmosphere 
models, respectively. The dashed line in the top panel is fixed to show a slope of unity. 
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Fig. 18. — Similar to Figure El after removing systematic differences by bringing the EPM 
distances to the SCM scale (all in Mpc). Only random differences between EPM and SCM 
distances remain. For the E96 case the discrepancies have a spread of 16%, while for the 
D05 case the scatter is 13%. The dashed line in the top panel is fixed to show a slope of 
unity. 
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Table 1: Telescopes and Instruments 



Telescope 


Instrument 


Phot/Spec 


CTIO 0.9 m 


CCD 


P 


YALO 1.0 m 


ANDICAM 


P 


YALO 1.0 m 


2DF 


S 


CTIO 1.5 m 


CCD 


P 


CTIO 1.5 m 


CSPEC 


s 


Blanco 4.0 m 


CSPEC 


s 


Blanco 4.0 m 


2DF 


s 


Blanco 4.0 m 


CCD 


p 


Swope 1.0 m 


CCD 


p 


du Pont 2.5 m 


WFCCD 


P/S 

/ 


du Pont 2.5 m 


MODSPEC 


s 


du Pont 2.5 m 


2DF 


s 


du Pont 2.5 m 


CCD 


p 


Baade 6.5 m 


LDSS2 


P/S 


Baade 6.5 m 


B&C 


s 


Clay 6.5 m 


LDSS2 


P/S 


ESO 1.5 m 


IDS 


s 


Danish 1.5 m 


DFOSC 


P/S 


MPI/ESO 2.2 m 


EFOSC2 


s 


NTT 3.6 m 


EMMI 


s 


ESO 3.6 m 


EFOSC 


s 


Kuiper 1.5 m 


CCD 


p 


Bok 2.2 m 


B&C 


s 


Lick 3.0 m 


Kast 


s 


Whipple 1.5 m 


FAST 


s 



Note. — Whether the instrument was used for photometry (P), spectroscopy (S), or both (P/S) is listed 
in Column 3. 



-51 - 



Table 2: Supernova Sample 



SN Name 


Host Galaxy 


a 


(J2000) 


<5(J2000) 






■^-host 


(s) b 


E(B - V) Gal c 


References 






[h 


m s 


] 


[° 


' 


" ] 










[mag] 




1991al 


LEDA 140858 


19 


42 


24 


00 


-55 


06 


23 








01525 


HP02 


0.051 


1 


1992af 


ESO 340-G038 


20 


30 


40 


20 


-42 


18 


35 








01847 


NED 


0.052 


1 


1992am 


MCG-01-04-039 


01 


25 


02 


70 


-04 


39 


01 








04773 


NED 


0.049 


1 


1992ba 


NGC 2082 


05 


41 


47 


10 


-64 


18 


01 








00395 


NED 


0.058 


1 


1993A 


anonymous 


07 


39 


17 


30 


-62 


03 


14 








02800 


NED 


0.173 


1 


1999br 


NGC 4900 


13 


00 


41 


80 


+02 


29 


46 








00320 


NED 


0.024 


2 


1999ca 


NGC 3120 


10 


05 


22 


90 


-34 


12 


41 








00931 


NED 


0.109 


2 


1999cr 


ESO 576-G034 


13 


20 


18 


30 


-20 


08 


50 








02020 


NED 


0.098 


2 


1999em 


NGC 1637 


04 


41 


27 


04 


-02 


51 


45 


2 





00267 


NED 


0.040 


2 


1999gi 


NGC 3184 


10 


18 


17 


00 


+41 


25 


28 








00198 


NED 


0.017 


3 


0210 


MCG +00-03-054 


01 


01 


16 


80 


-01 


05 


52 








05140 


NED 


0.036 


4 


2002fa 


NEAT J205221.51 


20 


52 


21 


80 


+02 


08 


42 








06000 


NED 


0.099 


1 


2002gw 


NGC 922 


02 


25 


02 


97 


-24 


47 


50 


6 





01028 


NED 


0.020 


4 


2002hj 


NPM1G +04.0097 


02 


58 


09 


30 


+04 


41 


04 








02360 


NED 


0.115 


1 


2002hx 


PGC 23727 


08 


27 


39 


43 


-14 


47 


15 


7 





03099 


NED 


0.054 


1 


2003B 


NGC 1097 


02 


46 


13 


78 


-30 


13 


45 


1 





00424 


NED 


0.027 


4 


2003E 


MCG-4-12-004 


04 


39 


10 


88 


-24 


10 


36 


5 





01490 


J09 


0.048 


1 


2003T 


UGC 4864 


09 


11 


11 


06 


+16 


44 


48 








02791 


NED 


0.031 


1 


2003M 


NGC 5374 


13 


57 


30 


65 


+06 


05 


36 


4 





01459 


J09 


0.027 


1 


2003bn 


2MASX J10023529 


10 


02 


35 


51 


-21 


10 


54 


5 





01277 


NED 


0.065 


1 


2003ci 


UGC 6212 


11 


10 


23 


83 


+04 


49 


35 


9 





03037 


NED 


0.060 


4 


2003cn 


IC 849 


13 


07 


37 


05 


-00 


56 


49 


9 





01811 


J09 


0.021 


1 


2003cx 


NEAT J135706.53 


13 


57 


06 


46 


-17 


02 


22 


6 





03700 


NED 


0.094 


4 


2003ef 


NGC 4708 


12 


49 


42 


25 


-11 


05 


29 


5 





01480 


J09 


0.046 


4 


2003fb 


UGC 11522 


20 


11 


50 


33 


+05 


45 


37 


6 





01754 


J09 


0.183 


4 


2003gd 


M74 


01 


36 


42 


65 


+15 


44 


20 


9 





00219 


NED 


0.069 


4, 5, 6 


2003hd 


MCG-04-05-010 


01 


49 


46 


31 


-21 


54 


37 


8 





03950 


NED 


0.013 


4 


2003hg 


NGC 7771 


23 


51 


24 


13 


+20 


06 


38 


3 





01427 


NED 


0.074 


1 


2003hk 


NGC 1085 


02 


46 


25 


76 


+03 


36 


32 


2 





02265 


NED 


0.037 


4 


2003hl 


NGC 772 


01 


59 


21 


28 


+19 


00 


14 


5 





00825 


NED 


0.073 


4 


2003hn 


NGC 1448 


03 


11 


36 


10 


-44 


37 


49 








00390 


NED 


0.014 


1 


2003ho 


ESO 235-G58 


21 


06 


30 


56 


-48 


07 


29 


9 





01438 


NED 


0.039 


1 


2003ip 


UGC 327 


00 


33 


15 


40 


+07 


54 


18 








01801 


NED 


0.066 


4 


2003iq 


NGC 772 


01 


59 


19 


96 


+18 


59 


42 


1 





00825 


NED 


0.073 


1 


2004dj 


NGC 2403 


07 


37 


17 


00 


+65 


35 


58 


1 





00044 


NED 


0.040 


7 


2004et 


NGC 6946 


20 


35 


25 


30 


+60 


07 


18 








00016 


NED 


0.342 


8 


2005cs 


NGC 5194 


13 


29 


53 


40 


+47 


10 


28 








00154 


NED 


0.035 


9, 10 



a Heliocentric host-galaxy redshift. 
b Source of host-galaxy redshift: HP02 
NASA/IPAC Extragalactic Database. 



Hamuv k Pmtol |2002h ; J09 = bones et all |2009l) ; NED 



: Galactic reddenings were obtained from ISchlegel et ah ( 1998h 



References. — (1) Calan/Tololo Super nova Program; (2) SO IRS; (3 )ILeonard et al.l (I2002bl) ; (4) Carnegie 
Type II S upernova Survey ( CAT S): (5) IVan Dvk et al. h00$l: (6) iHendrv et all J2005h : (7) Ivinko et al. 
d2006h : C8l lSahu et all (12009) : (9) IPastorello et al.l feoOflV floV lTsvetkov et al.l t00^ ~ 
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Table 3: Parameters from the Analytical Fitting to the V Light Curve 



SN Name 




a 




w 


Po 


too 


P 


Q 


R 


xl 




[mag] 


[JD-2, 448, 0001 

L ' ' J 


L 


lay] 


[mag/dayl 

L 0/ J J 


[mag] 


[magi 

L OJ 


[day] 


fdavl 

L 'J J 




1991al 


1 


.970 


519±3 


4 


.713 


0.009 


19. 


.061 


-0.244 


492.5 


28.44 


1.7 


1992af 


2 


.119 


861.1 


6 


.080 


0.024 


18. 


.094 


-2.614 


784.5 


90.07 


1.7 


1992am 


1 


.968 


947±5 


9 


.898 


0.006 


21 


.466 


0.296 


836.3 


37.61 


2.2 


1992ba 


1 


.976 


1014.0 


6 


.358 


0.008 


18. 


.033 


-0.081 


921.8 


16.35 


4.1 


1993A 


3 


.923 


1107±20 





.108 


0.009 


24. 


.229 


-0.392 


1009.5 


87.23 


9.8 


1999br 


1 


.937 


3396±20 


1 


.902 


0.013 


19. 


.901 


-2.319 


206.4 


162.85 


9.6 


1999ca 


2 


.841 


3373.2 


6 


.227 


0.011 


19. 


.975 


1.860 


268.4 


31.70 


1.8 


1999cr 


2 


.014 


3350±5 


7 


.939 


0.008 


20. 


.684 


-0.395 


294.0 


39.02 


2.6 


1999em 


1 


.978 


3590.1 


4 


.522 


0.009 


16. 


.297 


-0.403 


475.3 


64.56 


5.3 


1999gi 


2 


.131 


3645.4 


4 


.891 


0.006 


17. 


.470 


-1.941 


514.2 


6.96 


2.1 


0210 


2 


.120 


4592±5 


5 


.740 


0.009 


22. 


.768 


-0.601 


603.7 


62.48 


1.6 


2002fa 


1 


.701 


4577±10 


6 


.440 


0.010 


21. 


.624 


-0.906 


547.8 


54.02 


3.1 


2002gw 


2 


.131 


4661.2 


6 


.151 


0.008 


19. 


.691 


-0.451 


570.0 


57.69 


2.3 


2002hj 


1 


.968 


4660.3 


5 


.973 


0.014 


21. 


.470 


0.149 


570.3 


46.50 


4.8 


2002hx 


1 


.704 


4658.9 


5 


.884 


0.011 


20 


.802 


-0.532 


627.2 


47.42 


5.5 


2003B 


2 


.130 


4713.8 


6 


.172 


0.011 


18. 


.121 


-0.563 


641.4 


63.92 


1.9 


2003E 


2 


.677 


4766±5 


7 


.090 


0.009 


21. 


.563 


-0.542 


650.4 


42.14 


10.0 


2003T 


1 


.871 


4758.7 


3 


.671 


0.010 


21 


.534 


-0.176 


676.5 


46.69 


2.0 


2003bl 


3 


.886 


4805±5 


3 


.803 


0.011 


21. 


.893 


-1.918 


682.2 


162.52 


3.7 


2003bn 


2 


.237 


4813.7 


5 


.660 


0.007 


20. 


.076 


-0.168 


726.7 


25.60 


3.2 


2003ci 


3 


.105 


4818±5 


5 


.204 


0.012 


22. 


.687 


-0.705 


783.0 


59.46 


2.1 


2003cn 


3 


.524 


4804±5 


3 


.770 


0.019 


22. 


.587 


-0.654 


689.7 


129.94 


6.2 


2003cx 


1 


.980 


4828.5 


8 


.452 


0.006 


21. 


.519 


-0.492 


770.8 


67.70 


4.6 


2003ef 


2 


.724 


4870±5 


5 


.607 


0.009 


21. 


.157 


0.120 


847.4 


49.27 


4.6 


2003fb 


I 


885 


4874.4 


6 


.755 


0.010 


21 


.482 


-0.572 


811.1 


50.70 


3.6 


2003gd 


2 


.439 


4840.9 


2 


.210 


0.010 


17. 


.308 


1.245 


766.1 


42.03 


10.6 


2003hd 


2 


.105 


4952.9 


1 


.475 


0.012 


21. 


.767 


-0.207 


892.5 


53.77 


5.9 


2003hg 


2 


.684 


4999±5 


7 


.800 


0.009 


21. 


.259 


-0.199 


914.7 


19.66 


3.2 


2003hk 


3 


.097 


4961±3 


6 


.296 


0.011 


21. 


.591 


-0.180 


930.6 


21.75 


6.3 


2003hl 


3 


.376 


5005±5 


7 


.809 


0.008 


20. 


.973 


-0.103 


892.3 


39.14 


9.3 


2003hn 


2 


.010 


4963.7 


3 


.054 


0.013 


17. 


.466 


-0.231 


923.4 


52.79 


8.8 


2003ho 


2 


.319 


4920.9 


4 


.101 


0.009 


21. 


.703 


-0.159 


906.2 


4.43 


0.2 


2003ip 


2 


.676 


4946±20 


6 


.029 


0.015 


20. 


.321 


-0.531 


978.4 


71.97 


8.6 


2003iq 


1 


.604 


5020±5 


8 


.690 


0.010 


18 


.116 


-0.326 


937.8 


50.85 


2.3 


2004dj 


2 


.361 


5283±3 


5 


.134 


0.008 


14. 


.703 


-0.323 


174.7 


54.01 


0.9 


2004et 


1 


.850 


5394±5 


7 


.439 


-0.004 


15. 


.810 


1.969 


237.5 


145.32 


8.6 


2005cs 


2 


.883 


5666±10 


4 


.830 


0.004 


17. 


.822 


-0.191 


575.9 


10.03 


2.2 



* The minimum uncertainty of tpx is set to a conservative value of 2 days, and assigned to those SNe for which 
an uncertainty is not quoted. For three cases ipx is manually determined with an uncertainty of 20 days. 
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Table 4: Host- Galaxy Extinctions 



SN Name 


A v (spec) a 


Subclass 


A v (NaI D) b 


A V (V- 


-If 


1991al 


0.31(16) 


silver 


U.ol 


v UO ) 


-0.17 


[21) 


1992af 


1.24(31) 


coal 


n 1 7 


'1 ^\ 


-0.37 


[21) 


1992am 






n on 




0.52 


;23) 


1992ba 


0.43(16) 


silver 


n on 
u.uu 


'n^ 

v UO ) 


0.30 


;2i) 


1993A 


0.00(31) 


bronze 


n on 

u.uu 




0.06 


;25) 


1999br 


0.25(16) 


silver 


n nn 
u.uu 


'(\A \ 
v U4 ) 


0.94 


;25) 


1999ca 


0.12(31) 


coal 




v UO ) 


0.25 


;2i) 


1999cr 


0.47(31) 


coal 


n rq 


'91 


0.12 


;2i) 


1999em 


0.31(16) 


gold 


1 01 


'n^ 


0.24 


;2i) 


1999gi 


0.56(16) 


silver 


n ^n 
u.ou 


v UO J 


1.02 


;2i) 


0210 


0.31(31) 


bronze 


n nn 
u.uu 


V / 


0.31 


;36) 


2002fa 






n nn 
u.uu 


AV 


-0.35 


;26) 


2002gw 


0.40(19) 


silver 


n nn 

u.uu 


'n9"i 

\ ) 


0.18 


;22) 


2002hj 


0.16(31) 


bronze 


n nn 
u.uu 


v UO J 


0.24 


;22) 


2002hx 


0.16(25) 


coal 


n nn 
u.uu 


'1 R\ 


0.38 


;22) 


2003B* 


0.00(25) 


silver 


n 1 9 

U.1Z 


'n^ 

v UO ) 


-0.09 


'21) 


2003E 


1.09(31) 


coal 


n 71 

U. 1 1 


V<) 


0.78 


;23) 


2003T 


0.53(31) 


coal 


n 1 o 
u. iy 


'1 OA 


0.35 


'21) 


2003bl* 


0.00(16) 


gold 


n 1 1 

U. 1 1 


'1 n"; 


0.26 


[21) 


2003bn* 


0.09(16) 


silver 


n nn 
u.uu 


'n^ 

v UO J 


-0.04 


[21) 


2003ci 


0.43(31) 


coal 


n nn 
u.uu 


v ZO ) 


0.78 


;23) 


2003cn* 


0.00(25) 


gold 


n nn 

u.uu 


'noA 


-0.04 


'23) 


2003cx 


0.65(25) 


coal 






-0.27 


;25) 


2003ef 


1.24(25) 


gold 


i ztn 


'i 9"; 


0.98 


;2i) 


2003fb 


0.37(31) 


coal 


n ka 




1.24 


;23) 


2003gd 


0.40(31) 


coal 


0.00 


'04^1 
v ui; 


0.33 


;2i) 


2003hd 


0.90(31) 


coal 


0.74 


;27) 


0.01 


;2i) 


2003hg 






2.29 


;i2) 


1.97 


;24) 


2003hk 


0.65(31) 


coal 


1.69 


;20) 


0.44 


;25) 


2003hl 


1.24(25) 


gold 


1.84 


;o9) 


1.72 


;23) 


2003hn 


0.59(25) 


coal 


0.64 


;o8) 


0.46 


;2i) 


2003ho 


1.24(31) 


bronze 


1.28 


;io) 


2.19 


;2i) 


2003ip 


0.40(31) 


coal 


0.42 


'08) 


0.56 


;22) 


2003iq 


0.37(16) 


silver 


0.91 


;o4) 


0.25 


;22) 


2004dj 


0.50(25) 


coal 


0.26 


;o6) 


-0.09 


;23) 


2004et 


0.00(25) 


coal 


1.17 


;o2) 


0.13 


;27) 


2005cs 










0.72 


;30) 



a Spectroscopic reddenings assuming Ry = 3.1. Additionally, the third column lists the subclass defined from 
the criteria given in § 14.11 

6 Equivalent width of Na I D measured by us and converted to Ay with the relationship of Barbon ct al 
(|l990h using R v = 3.1. 
c This paper. 

* SNe selected to determine the intrinsic color. 
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Table 5: Magnitudes, Expansion Velocities, and V — I Colors for Day -30 



SN Name 






B 






V 




I 


^Fe II 


V - 


-I 




±187 [km s- 1 ] 




[mag] 




[mag] 




[mag] 


[km s- 1 ] 


[mag] 


1991al 


4480 


17 


944 


'036"! 


16 


917(029) 


16 


333(020) 


5328(202) 


589 




1992af 


5359 


18 


024 


'038) 


17 


096(029) 


16 


706(053) 


4529(212) 


0.511 


'027) 


1992am 


14007 


20 


443 


'1 1 7) 


19 


053 038) 


18 


218 033) 




830 


'046"! 

. U4:U J 


1992ba 


1245 


16 






15 


657(031) 


14 


886(022) 


2237(068) 


U. M U 


'091 ^ 

V / 


1993A 


8907 


20 


uou 


'1 97\ 


19 


309(089) 


18 


657(057) 




fi78 

U.UI o 




1999br 


1285 


19 


010 


'022) 


17 


582(013) 


16 


574(008) 


1127(205) 


1.028 


'011) 


1999ca 


3108 


17 


901 


'071) 


16 


489(055) 


15 


735 034) 




0.755 


'025) 


1999cr 


6363 


19 


435 

I'.) O 


'067") 

, UU 1 > 


18 


418(051) 


17 


723(035) 


3095(207) 


703 

U . i UO 




1999em 


670 


15 


331 


'034) 


13 


998(023) 


13 


245(017) 


2727(056) 


0.752 


'019) 


1999gi 


831 


16 


554 


'032) 


15 


060(019) 


14 


011(022) 


2725(061) 


1.058 


'025) 


0210 


15082 


21 


955 


'044) 


20 


572(028) 


19 


733(035) 


4923(206) 


0.780 


'110) 


2002fa 


17847 


21 


222 


'0721 


20 


243(052) 


19 


709(088) 


4176(061) 


518 

U . O J- o 


'065"i 


2002gw 


2878 


18 


575 

U 1 O 


'090") 

UZ(U J 


17 


491(029) 


16 


752(014) 


2669(063) 


726 


'035"l 

. UOU J 


2002hj 


6869 


19 


705 

1 UU 


'084 s ! 


18 


582(061) 


17 


939(042) 


3657(065) 


751 


'031 "l 

V / 


2002hx 


9573 


20 


^98 




19 


164(035) 


18 


360(028) 


3960(203) 


u . ouu 


V / 


2003B 


1105 


17 


078 


'022) 


15 


966(018) 


15 


341(015) 


2795(067) 


0.620 


'017) 


2003E 


4380 


19 


659 


'046) 


18 


368(026) 


17 


472(024) 


2869(203) 


0.965 


'024) 


2003T 


8662 


20 


541 


'098") 


19 


227(020) 


18 


423(022) 


2803(203) 


7Q5 

U . i UO 


'094"l 

V / 


2003bl 


4652 


20 


J.OO 


'049^ 


18 


961(020) 


18 


230(019) 


1475(202) 


U. ( oo 


'09^ 

. \J£0 I 


2003bn 


4173 


18 


KIR 
uoo 




17 


403(034) 


16 


769(024) 


2719(202) 


U . U4: 1 




2003ci 


9468 








19 


634(035) 


18 


648(029) 


2876(077) 


0.963 


'034) 


2003cn 


5753 


19 


852 


'024) 


18 


827(015) 


18 


182 014) 


2510(207) 


0.640 


'017) 


2003cx 


11282 


20 


41 7 


'1 noi 

v / 


19 


513(050) 


19 


050(067) 


3734(205) 


551 

U . OO -L 


'058"l 

, UOO 1 


2003ef 


4504 


19 


304 


'032 s ! 


17 


837(024) 


16 


792(017) 




1 043 

_L •\J L ±0 


'01 5") 

V / 


2003fb 


4996 


20 


521 


'122) 


19 


085(092) 


17 


940(056) 


3120(206) 


1.148 


'040) 


2003gd 


359 


15 


208 


;045) 


13 


965(036) 


13 


167(023) 


2976(230) 


0.786 


'019) 


2003hd 


11595 


20 


462 


'031) 


19 


320(014) 


18 


658(018) 


3741(069) 


0.661 


;022) 


2003hg 


3921 


20 


247 


'091) 


18 


067(040) 


16 


603(025) 


3398(211) 


1.435 


'027) 


2003hk 


6568 


19 


478 


;036) 


18 


201(024) 


17 


380(023) 


3618(087) 


0.828 


^026) 


2003hl 


2198 


19 


124 


;062) 


17 


219(038) 


15 


806(027) 


2354(063) 


1.336 


'027) 


2003hn 


1102 


16 


416 


;020) 


15 


153(011) 


14 


307(015) 


3121(074) 


0.836 


'018) 


2003ho 


4134 


20 


779 


;084) 


18 


946(026) 


17 


419(016) 


4152(080) 


1.523 


^022) 


2003ip 


5050 


18 


894 


;046) 


17 


549(035) 


16 


664(025) 


3813(204) 


0.876 


'021) 


2003iq 


2198 


17 


396 


;049) 


16 


124(037) 


15 


362(024) 


3482(063) 


0.756 


'018) 


2004dj 


180 


12 


973 


;074) 


12 


046(035) 


11 


416(033) 


2725(202) 


0.621 


'043) 


2004et 


-133 


13 


683 


^220) 


12 


099(171) 


11 


378(104) 


2901(202) 


0.707 


'070) 


2005cs 


635 


16 


099 


;062) 


14 


749(051) 


13 


824(055) 




0.957 


^067) 



* The Velocity Calculator tool at the NASA/IPAC Extragalactic Database webpage computes the CMB red- 
shift (czcmb in km s _1 ) given the coordinates and the heliocentric redshift of an object. Its error corresponds 
to the uncertainty in the determination of the Local Group velocity (187 km s _1 ). 

Note. - All of the data in this table have been interpolated to day —30 with errors accounting for 
the interpolation, the error in tpT, and the instrumental error. BVI corresponds to apparent magnitudes 
corrected for Galactic absorption and if-terms, but uncorrected for host-galaxy extinction. The V — I colors 
do not exactly match the V minus / magnitude difference, since light and color curves are interpolated 
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Table 6: Fitting Parameters from the Hubble Diagrams 



Reddening 










estimator 


filter 


a 




zp 


(V -iy 


B 


3.27 ±0.37 




-0.24 ±0.08 




V 


3.01 ±0.36 




-1.36 ±0.08 




I 


2.99 ±0.36 




-2.02 ±0.08 


Spectral 


B 


4.77 ±0.47 




-0.60 ±0.11 


analysis 13 


V 


4.22 ± 0.45 




-1.60 ±0.10 




I 


3.67 ±0.43 




-2.11 ±0.10 


(^-■0 + 


B 


3.50 ±0.30 


2.67 ±0.13 


-1.99 ±0.11 


variable Ry c 


V 


3.08 ±0.25 


1.67 ±0.10 


-2.38 ±0.09 




I 


2.62 ±0.21 


0.60 ±0.09 


-2.23 ±0.07 



a Reddenings from § 13.21 

6 Spectroscopic reddenings. 

c Leaving R v as a free parameter. 

Note. — These are the results of fitting Equations (TT21 and (|T5j) to the data in Tables 2] and [5] 
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Table 7: H Q Calculations 





M + a log(t;/5000) - /3{V -■ 


I) [mag] 


Ho 


[km s 1 Mpc 


- 1 ] 


SN 


B V 


I 


H (B) 


H (V) 


Ho(I) 


1999em 


-17.94 ±0.20 -18.41 ±0.20 


-18.24 ±0.19 


64.7 ±6.8 


62.2±6.1 


62.9±6.0 


2004dj 


-17.08 ±0.30 -17.28 ±0.27 


-17.13 ±0.26 


95.7±14.0 


104.5 ±13.7 


104.8 ±12.9 


Average 


-17.67 ± 0.39 -18.02 ±0.53 


-17.84 ±0.53 


71 ±12 


69 ±16 


70 ±16 



Note. — The uncertainties in (Ho) and (M corr ) (the corrected absolute magnitude for BVI in the second 
column) correspond to the weighted dispersion in each pair of values, which are a factor of 2-3 greater than 
the formal errors. This is a conservative uncertainty which takes into account the larger than expected 
discrepancy in absolute magnitude between the two objects. 
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Table 8: Distance Moduli [mag] 
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U4^(J4 J 


34 

04 
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00 


31 (^4"! 
01 ^04 J 
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00 


98(^3"! 


35 


30(38) 


9003rn 
ZUUOC11 


TP 8zLQ 
ivy o^ty 


34 
04 


77(49"! 

V / 


34 

04 


8fi( c ; c ;"\ 


34 

04 


8^(^4"l 


34 


81(28) 


9003rv 


NFAT T13^7flfi ^3 


3fi 
ou 


1 7(44"! 

1 / \^44 J 


3fi 
ou 


93( c ; c ;"! 

zo yoo j 


3fi 
ou 


93(^4"! 

ZO\^04 J 


36 


20(29) 


9003fh 

ZUUOIU 


TTPP 1 1 ^99 
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34 

04 


41 (44 "l 


34 

04 


OU^OO J 


34 

04 


^(^4"! 


34 


49(29) 
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35 

OO 
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35 

OO 


85(54"! 


35 

OO 


77(53") 

1 I I OO J 


35 


86(28) 


2003hg 


NGC 7771 


33 


50(42) 


33 


18(54) 


33 


14(54) 


33 


31(28) 


2003hk 


NGC 1085 


34 


45(40) 


34 


41(54) 


34 


35(53) 


34 


41(28) 


2003hl 


NGC 772 


32 


08(41) 


32 


01(54) 


32 


00(53) 


32 


04(28) 


2003hn 


NGC 1448 


31 


14(40) 


31 


15(54) 


31 


11(53) 


31 


13(27) 


2003ho 


ESO 235-G58 


34 


10(41) 


34 


18(54) 


34 


13(53) 


34 


13(28) 


2003ip 


UGC 327 


33 


81(41) 


33 


75(54) 


33 


67(54) 


33 


76(28) 


2003iq 


NGC 772 


32 


50(40) 


32 


40(54) 


32 


34(53) 


32 


43(28) 


2004dj 


NGC 2403 


28 


06(43) 


28 


22(55) 


28 


19(54) 


28 


14(29) 


2004et 


NGC 6946 


28 


64(50) 


28 


21(58) 


28 


17(55) 


28 


37(31) 



Note. — These values were calculated according to the approach discussed in § 14.51 The mean value is a 
weighted average of the individual values for each filter. 



Table 9: EPM and SCM Distances 



SN Name 


<^E96 a 


0?D05 b 


dscM 




[Mpc] 


[Mpc] 


[Mpc] 


1992ba 


16.4 (2.5) 


27.2 


(6.5) 


18.5 (2.4) 


1999br 




39.5 


(13.5) 


23.6 (3.4) 


1999em 


9.3 (0.5) 


13.9 


(1.4) 


10.0 (1.3) 


1999gi 


11.7 (0.8) 


17.4 


(2.3) 


12.6 (0.6) 


2002gw 


37.4 (4.9) 


63.9 


(17.0) 


48.1 (6.2) 


2003T 


87.8 (13.5) 


147.3 


(35.7) 


108.1 (14.0) 


2003bl 




92.4 


(14.2) 


65.2 (9.0) 


2003bn 


50.2 (7.0) 


87.2 


(28.0) 


52.5 (6.8) 


2003hl 


17.7 (2.1) 


30.3 


(6.3) 


25.6 (3.3) 


2003hn 


16.9 (2.2) 


26.3 


(7.1) 


16.8 (2.1) 


2003iq 


36.0 (5.6) 


53.3 


(17.1) 


30.6 (4.0) 



EPM distances from atmosphere models by Eastman et al. ( 19961) : see Jones et al. ( 2009 ) 



EPM distances from atmosphere models by Dessart fc Hillier ( 20051) : see Jones et al. ( 2009 ). 
Note. — Uncertainties are given in parentheses. 



